Creep of current-driven domain-wall lines: intrinsic versus extrinsic pinning 
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We present a model for current- driven motion of a magnetic domain-wall line, in which the dynam- 
ics of the domain wall is equivalent to that of an overdamped vortex line in an anisotropic pinning 
potential. This potential has both extrinsic contributions due to, e.g., sample inhomogeneities, and 
an intrinsic contribution due to magnetic anisotropy. We obtain results for the domain-wall velocity 
as a function of current for various regimes of pinning. In particular, we find that the exponent char- 
acterizing the creep regime depends strongly on the presence of a dissipative spin transfer torque. 
We discuss our results in the light of recent experiments on current-driven domain-wall creep in 
ferromagnetic semiconductors, and suggest further experiments to corroborate our model. 

PACS numbers: 72.25.Pn, 72.15.Gd, 75.60.Ch, 85.75.-d 
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I. INTRODUCTION 

The driven motion of line defects through a dis- 
ordered potential landscape has attracted considerable 
attention, for example in the context of vortices in 
superconductors,^ wetting phcnomenaf^, crack fronts,'^ 
and domain walls in ferromagnets.'*-i^ The competition 
and interplay among the elasticity of the line, the pin- 
ning forces due to the disorder potential, and thermal 
fluctuations, lead to a wealth of physical phenomena. 
Topics discussed are, for example, the universality class 
of the roughening of the line, the nature of the pinning- 
dcpinning transition at zero temperature,^ and the slide, 
depinning, and creep regimes of motion of the line that 
occur for decreasing driving fieldiiii 

The creep regime has been observed experimen- 
tally with field-driven motion of domain walls in 
ferromagnets4:^ This low-field regime is characterized by 
a nonlinear dependence of the domain-wall drift velocity 
{X) on the external magnetic field -ffoxt, given by 



{X) (X exp < — 



Er. 



(1) 



where Ec is a characteristic energy scale, and He a crit- 
ical field. The thermal energy is denoted by ksT and 
the exponent = (2C — l)/{2 — Q is given in terms 
of the equilibrium wandering exponent C of the static 
linejii^ii The phenomenological creep formula [Eq. ([T])] 
is an Arrhenius law in which the energy barrier diverges 
for vanishing driving field, i.e., is "glassy". The underly- 
ing assumption is that there is a characteristic length 
scale that determines the displacement of the domain 
wall line. The validity of Eq. ((T|) has been confirmed 
both numerically^ and with functional renormalization 
group methods!^ It turns out that Eq. |T]) is also valid in 
situations where roughening plays no role. For example, 
for a d-dimcnsional manifold driven through a periodic 
potential in d -I- 1 dimensions we have /i/ ~ d — 1 (for 
d > 2)fi Moreover, in the regime where the line defect 



moves via variable-range hopping, we have that /i/ — 1/3 
if the motion is in two dimensions iii^iii 

In addition to magnetic-field driven motion, a lot of 
recent theoretical and experimental research has been 
devoted to manipulating domain walls with electric 
curreniii^ii^iiiii^ii^iiiii^ii^i^i2i via so-called spin transfer 
torques i^i^i^iiS^ Domain wall motion driven by a cur- 
rent is quite distinct from the field-driven case. For ex- 
ample, it has been theoretically predicted that, in certain 
regimes of parameters, the domain wall is intrinsically 
pinned at zero temperature, meaning that there exists a 
nonzero critical current even in the absence of disorder 
In clean samples, the phenomenology of current-driven 
domain wall motion turns out to crucially depend on the 
ratio of the dissipative spin transfer torque parameter— 
13 and the Gilbert damping constant Although 
theoretical prediction a^^'^^i^^'^" indicate that, at least 
for model systems, this ratio differs from one, it turns 
out to be difficult to extract its precise value from ex- 
periments on current-driven domain wall motion, to a 
large extent because disorder and nonzero-temperature 
effects^i^ complicate theoretical calculations of the do- 
main wall drift velocity for a given current. This is the 
first motivation for the work presented in this paper. 

Previous work on current-driven domain wall motion 
at nonzero temperature focused on rigid domain walls. 
Tatara et al^ found that ln(X) was proportional to 
the current density j. Discrepancy of this result with 
experimentSf^L that did not observe linear dependence 
of wall velocity on current, motivated the more system- 
atic inclusion of nonzero-temperature effects on rigid- 
domain wall motion by Duine et al.,^^ who found that 
ln(X) cx ^/j in certain regimes. Although the latter 
was an important step in qualitatively understanding 
the experimental results of Yamanouchi et al.,^^'^^ a de- 
tailed understanding of these experiments is still lack- 
ing and this is the second motivation of this paper. For 
completeness, we mention also the theoretical work by 
Martinez et alM--^ who considered thermally-assisted 
current-driven rigid domain wall motion in the regime of 
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Fig. 1: Duine/Morais Smith 



FIG. 1: Mapping of current-driven domain wall dynamics to 
that of a vortex line. The position of the domain wall X{z, t) 
and its chirality <j)o{z,t) become the position {ux,Uy) of the 
vortex via {ux,Uy) = {X/X, 4>o)- The potential landscape for 
this vortex is in general anisotropic. In particular, the tilting 
in the Ua;-direction is set by the external magnetic field and 
the dissipative spin transfer torque. The tilting in the Uy- 
direction is determined by the reactive spin transfer torque. 



large anisotropy, where the chirality of the domain wall 
plays no role and the pinning is essentially dominated by 
extrinsic effects. Furthermore, Ravelosona et al^ have 
observed thermally-assisted domain wall depinning, and 
Laufenberg et al.^'' have determined the temperature de- 
pendence of the critical current for depinning the domain 
wall. 

In this paper we present a model for a current-driven 
elastic domain-wall line moving in one dimension in the 
presence of disorder and thermal fluctuations. A cru- 
cial ingredient in the description of current-driven mo- 
tion is the chirality of the domain wall, which acts like 
an extra degree of freedom. This enables a reformu- 
lation of current-driven domain wall motion as a vor- 
tex line moving in an anisotropic potential in two di- 
mensions (see Fig. [T]), which we present in detail in 
Secini Using this physical picture, we analyze in Sec. lIIII 
the different regimes of pinning within the framework 
of collective pinning theory,^ We present results on the 
velocity of the domain-wall line as a function of cur- 
rent, both in the regime where intrinsic pinning due 
to magnetic anisotropy dominates, and in the extrinsic- 
pinning-dominated regime. Finally, we discuss in Sec. II VI 
our theoretical results in relation to recent experiments 
on current-driven domain walls in GaMnAs."^'^ Although 
these experiments remain, in our opinion, not fully un- 
derstood, we suggest that they may be explained by as- 
suming a specific form of the pinning potential for the 
domain-wall line. We suggest futher experiments that 
could corroborate this suggestion. 



II. DOMAIN WALL AS A VORTEX LINE 

The equation of motion for the magnetization direction 
n in the presence of a transport current j is, to lowest 
order in temporal and spatial derivatives, given by 

(^^ + vs • V j f2 - O X (H + He,t + ry) = 

-aa^ X ( I- + • V ) f2 . (2) 

\ot ac J 

The left-hand side of this equation contains the reactive^ 
spin transfer torqu^ proportional to the velocity Vg = 
Pj/ (eps)- The latter velocity characterizes the efficiency 
of spin transfer. Here, P is the polarization of the current 
in the ferromagnet, e is the carrier charge, and the spin 
density is denoted by ps = with a the lattice con- 
stant. The other terms on the left-hand side of Eq. ^ 
describe precession around the external field Hoxt and 
the effective field H = —6E[i},]/{h6Cl), which is given 
by a functional derivative of the energy functional E[^l] 
with respect to the magnetization direction. The stochas- 
tic magnetic field tj incorporates thermal fluctuations, 
and it has zero mean and correlations determined by the 
fluctuation-dissipation theorem'^^ 

(r;,(x,<)r,,,(x',t')> = ^^^^,5(t - i')a'^(x - x')<5..' . 

(3) 

It can be shown that the above still holds in the pres- 
ence of current, at least to first order in the applied 
electric field'^*' that drives the transport current. The 
fluctuation-dissipation theorem also ensures that in equi- 
librium the probability distribution for the magnetization 
direction is given by the Boltzmann distribution VlTl] cx 
exp{-E[n]/kBT}. The right-hand side of Eq. © con- 
tains only dissipative terms. The Gilbert damping term 
is proportional to the damping parameter ac, and the 
dissipative^*' spin transfer torque is characterized by the 
dimensionless parameter /3fi^ 

We consider a ferromagnet with magnetization direc- 
tion rj = (sin 6 cos 4>, sin sin cj), cos 6) that depends only 
on the the x and z-direction. In addition, we take the 
current in the a;-direction and the external magnetic fleld 
in the z-direction. The size of the ferromagnetic fllm in 
the a direction is denoted by (a G {x,y, z}) and we 
assume that Ly <^ L^. The latter assumption allows us 
to model the domain wall as a line. Furthermore, we 
take the ferromagnet to have an easy z-axis and hard y- 
axis, with anisotropy constants K and Ki_, respectively. 
The spin stiffness is denoted by J. With these assump- 
tions, static domain walls have a width A = J / K , 
and are, for the simplest model to be discussed in more 
detail below [see Eq. (|10p ]. described by the solutions 
6'o(x) = cos^^ [tanh(x/A)] and 0(x) — 0. To arrive at 
a description of the dynamics of the domain wall, we 
use two collective coordinates which may depend on the 
z-coordinate, so that the domain wall is modelled as a 
line. The collective coordinates are the position of the 
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wall X{z,t) and the chirality (l)o{z,t). The latter deter- 
mines the sense in which the magnetization rotates u pon 
going through the domain wall. The result of Ref. [31[ 
is straightforwardly generalized to the case of a domain 
wall line. This amounts to solving Eq. ([2]) variation- 
ally with the ansatz 0dw(x, t) ~ 0q {{x — X{z,t))/X) and 
0dw(x, t) — (j)oiz, t), which yields the equations of motion 



at X dt 
1 dX 900 



-a^ SV (3vs „ , , 



6V 



where the domain- wall energy 

V[X,(bo] = E[0d^, 



WdwJ 



(4) 



(5) 



and the stochastic forces are determined from 



{r^^{z,t)r^4z',t')) = {r^x{z,t)rjx{z',t')) 



( 



\X^L, 



X 



5{t-t'). (6) 



The above equations are derived using a variational 
method for stochastic differential equations based on 
their path-integral formulation.-SliiS Their validity is con- 
firmed a posteriori by noting that in equilibrium the 
probability distribution function for the position and chi- 
rality of the domain wall is the Boltzmann distribution. 
That is, the Fokker-Planck equation for the probability 
distribution V [X, 0o] of the domain wall position and chi- 
rality that follows from Eqs. ^ and is given by^ 



(1 



2s 2hLyX^ dr[X,(f)a] 



dz 



agX 



dt 
SV 



SV 



+ 



SX{z) 
dz S 

S^ 



X- 



SX{z) SMz 
SV 



V[X,i 



(5 V 

"GTTT7T ]'P[X,(l>a] 



'S(j)o{z) 



+aGkBT / — 



A \S(t,l{z) SX^{z) 



V[X,M (7) 



Upon insertion of the Boltzmann distribution 
Vcq[X,(l)Q] oc exp{—V[X,(f)Q]/{kBT)} into this equa- 
tion, one straightforwardly verifies that it is indeed a 
stationary solution. 

Rewriting the equations of motion for the domain wall 
position and chirality in terms of the dimcnsionlcss coor- 
dinate u(z, t) = {X{z, t)/X, (j)o{z, t)), we find from Eq. (g]) 
that the domain wall is described by 



(z,t) = -aGUa{z,t) 



SV[\i] 

SUa{z, t) 



ilc{z,t) , (8) 



with e^a' the two-dimensional Levi-Civita symbol. 
(Summation over repeated indices a, a' € x, y is implied. 
Note that ry^, = rjx^tf, for a = x,y. ) The above equa- 
tion of motion [Eq. ([5])] corresponds to the overdamped 



limit of vortex-line dynamics in an anistropic potential 
V[vl\. The left-hand side of Eq. ^ corresponds to the 
Magnus force on the vortex. We emphasize that a mass 
term is missing, indicating we are indeed dealing with 
the overdamped limit of vortex motion. (Note that the 
mass of the fictitious vortex is not related to the Doring 
domain wall mass"*^ that arises from eliminating the chi- 
rality from the domain-wall description, which is valid 
provided the latter is small 4^ As the dynamics of the 
domain-wall chirality is essential for current-driven do- 
main wall motion, this latter approximation is not suffi- 
cient for our purposes.) The right-hand side of the equa- 
tion of motion contains a damping term proportional to 
ac, and a term representing thermal fluctuations. The 
force is determined by the potential 



V 



a^V[Xux,Uy 
2hLyX^ 



dz 

T 



A 



, . . (9) 

The tilting of this potential in the U3;-direction is deter- 
mined by the parameter (3, the current v^^ and the ex- 
ternal field i?oxt ■ The tilting in the Uj^-direction is deter- 
mined only by the current. The model in Eqs. ^ and ([5]), 
illustrated in Fig.[Tl is the central result of this paper. In 
the following section, we obtain results from this model 
for the domain wall velocity in different regimes of pin- 
ning, specializing to the case of current-driven domain 
wall motion (i/ext = 0). 



III. DOMAIN WALL CREEP 

In this section, we obtain results for the average drift 
velocity of the domain wall as a function of applied cur- 
rent. First, we discuss the situation without disorder, 
hereafter we incorporate the effects of disorder. 



A. Intrinsic pinning 

In this subsection, we make two assumptions that are 
not related. First, we consider a homogeneous system, 
i.e., a system without disorder potential l^in = 0. Sec- 
ond, we take 13 — 0. In this case, the domain wall is 
intrinsically pinned. This comes about as follows. For 
the magnetic nanowire model discussed in the previous 
section, the energy functional in the clean limit is given 
by 
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(10) 
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Upon insertion of the domain wall ansatz into the above 
energy functional, we find that 



J fdvL 



K 



X 



(11) 

the 



Because the above potential is independent of z 
domain wall remains straight at zero temperature, i.e., 
du/dz = 0. By solving the equations of motion in 
Eq. ([5]) for the potential in Eq. (fTT|) at zero temperature 
and for a straight domain wall, one finds that (|u|) oc 

■\J — {XK±/2Kf so that the domain wall is pinned up 
to a critical current given by jc = XK±eps/2hP. (The 
brackets (• • •) denote time and thermal average.) This in- 
trinsic pinning is entirely due to the anisotropy energyj^^ 
determined by K±, and does not occur for field-driven 
domain-wall motion, or current-driven domain-wall mo- 
tion with /3 7^ 0. Physically, it comes about because, for 
the model of a domain wall that we consider here, the 
reactive spin transfer torque causes the magnetization to 
rotate in the easy plane. This corresponds to an effective 
field that points along the hard axis. Because the Gilbert 
damping causes the magnetization to precess towards the 
effective field, the current tilts the magnetization out of 
the easy plane. This leads to a cost in anisotropy en- 
ergy which stops the drift motion of the domain wall if 
the current is too small. By solving the equations of mo- 
tion for the potential in Eq. (jlip at nonzero temperature 
in the limit of a straight wall, one recovers the result of 
Ref. [H. 

At nonzero temperature, the domain wall is no longer 
straight. Since only the chirality is important, our 
model for current-driven domain-wall motion in Eq. pip 
then corresponds to the problem of a string in a tilted- 
washboard potential, that has been studied before"'^ in 
different contexts. At nonzero temperature the string 
propagates through the tilted- washboard potential by nu- 
cleating a kink-antikink pair in the z-direction of the 
domain- wall chirality 4'o{z^ t). The kink and antikink are 
subsequently driven apart by the tilting of the potential, 
which results in the propagation of the string. 

In the limit when the current is close to the critical 
one, a typical energy barrier is determined by the com- 
petition between the elasticity of the string and the tilted 
potential^ For {jc — j)/ jc «C 1 the cosine in the energy 
functional in Eq. (fTTj) may be expanded around one of 
its minima, which yields 



d5i 




dz 



'-If 



(12) 



where we have omitted an irrelevant constant. In the 
above expression, 5uy denotes the displacement from the 
minimum. Note that we have dropped the dependence of 



the potential on Ux which is allowed because the potential 
is not tilted in the Ma;-direction (provided that (3 — 0). 

The potential in Eq. ((T2)) 
(by construction) and 



has a minimum for 5u^™ 



a maximum 



for 



{j I j c)'^ / Xh. The pinning potential energy 
barrier, i.e., the pinning potential evaluated at the maxi- 
mum, scales as AF cx [1 — (j/jc)^]^^^- Consider now the 
situation that a segment of length L of the string is dis- 
placed from the minimum and pinned by the maximum 
of the potential. The length L is then determined by the 
competition between the elastic energy ~ J{Su'^'^^/L)^, 
that tends to keep the domain wall straight, and the pin- 
ning potential Al^. Equating these contributions yields 
for the length L that 



L oc 



1 - 



-1/4 



(13) 



The typical energy barrier that thermal fluctuations have 
to overcome to propagate the domain wall is then given 
by evaluating Eq. (fT2|) for a segment of this length. This 
yields a typical energy barrier oc [1 — (j/jc)^]^^'*- Putting 
these results together and assuming an Arrhenius law, 
we find that the domain wall velocity is 



ln(|u|) cx - 



1 JLy 



K 



jc 



-1 5/4 



(14) 



for {jc-j)/jc < 1. 

In the regime of small currents j <JC jc, the typi- 
cal energy barrier depends only logarithmically on the 
current,*'* so that (|u|) oc j. This latter result for 
the domain-wall velocity is the same as found from 
a treatment of rigid domain-wall motion at nonzero 
temperature^* in the limit that j ^ jc- This is under- 
stood by noting that in the limit of vanishing current 
the elasticity of the domain-wall line does not enter the 
expression for the typical energy barrier,— the so-called 
"thin- wall" limit 4^ 



B. Extrinsic pinning 

We now add extrinsic pinning, i.e., a disorder potential 
Vpin to the potential in Eq. (fTTj) . Following Ref. [l2| 
we assume, in first instance, that it only couples to the 
position of the domain wall and not to its chirality Uy. 
This assumption is made mainly to simplify the problem. 
Considering now the general case that also /3 ^ 0, we have 
that 



2h 



du 

dz 



— — cos 2Uy 

Ah " 



+ ^^pin(Ma;, z) + (i^U^ + 



(15) 
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We estimate a typical energy barrier using the collective 
pinning theoryjiil Therefore, we assume that we are in 
the regime where the pinning energy grows sublinearly 
with the length of the wall, and that there exists a typ- 
ical length scale L at which domain- wall motion occurs. 
(Note that we take L dimensionless since the coordinate 
u is dimensionless.)^ The energy of a segment of this 
length that is displaced is given by 

2 

E{L) = eei^ + (3jLu.,. + jLuy . (16) 

The first term is the elastic energy with ed = J/2hX^. 
The second and third term correspond to the dissipa- 
tive and reactive spin transfer torques, respectively. Note 
that since the dissipative spin transfer torque acts like an 
external magnetic field we are able to incorporate it in the 
above energy. The potential Vp-miux, z) leads to roughen- 
ing in the direction. Following standard practicefii^'i 
we assume a scaling law Ux{L) = UxqL'', with C the equi- 
librium wandering exponent, already mentioned in the 
introduction, and u^o a constant. The displacement in 
the Uy direction is not roughened, because we have as- 
sumed that Vpiniuxjz) does not depend on u^, i.e., the 
domain-wall chirality. Rather, the displacement in this 
direction is determined by the minima of the potential in 
Eq. pip and we have that Uy = Uyo independent of L for 
j <C jc- Note that in this limit the elastic energy due to 
displacement in the Uy-direction can also be neglected;^ 
Hence, we find that 

E{L) = e^ml^L^^-^ + p^UxoL'^^' + jLuyo . (17) 

Minimizing this expression with respect to L then leads 
to a typical energy barrier. Assuming an Arrhenius 
lawfii^ii we find for the domain wall velocity 

For (3 = Owe have that fic = (2C- l)/(2 - 2C). For /? 7^ 
we find /ic = (2C— l)/(2 — C)- In particular, for C = 2/3, 
applicable to domain walls in ferromagnetic metals,^ we 
have /ic — 1/2 for /? = 0, and fic = 1/4 for (3^0. Since 
the dissipative spin transfer torque, proportional to /3, 
acts like an external magnetic field on the domain wall 
[see Eqs. ([5])], we recover the usual results for field-driven 
domain wall motion'* from our model. This result is also 
understood from the fact that an external magnetic field 
does not tilt the domain wall potential in the chirality di- 
rection, as opposed to a current, so that the domain wall 
chirality plays no role in field-driven domain wall creep. 
We note that Eq. (jj]), or, equivalently, Eq. ([5]), contains 
a description of Walker breakdown*^ in the clean zero- 
temperature limit and is also able to describe the transi- 



tion from the creep regime to the regime of precessional 
field-driven domain-wall motion observed recently.— 
IV. DISCUSSION AND CONCLUSIONS 

In very recent experiments on domain walls in the 
ferromagnetic semiconductor GaMnAs, Yamanouchi et 
al.^ have observed field-driven domain- wall creep with 
exponent ~ 1, and current-driven creep with — 
1/3, over 5 orders of magnitude of domain- wall veloc- 
ities. The fact that these two exponents are different 
could imply that /? is extremely small for this material. 
For /3 = and the specific pinning potential discussed 
in the previous section it is, however, impossible to find 
a single roughness exponent that yields both fij = 1 
and fic = 1/3. (Note that the theoretical arguments in 
Ref. ^3§\ give ^i/ = 1 and /ic = 1/2.) 

Although it is extremely hard to determine the mi- 
croscopic features of the pinning potential, we empha- 
size that if pinning is not provided mainly by point-like 
defects (as considered in this paper and argued by Ya- 
manouchi et al. to be the case in their experiments,'^'^) 
but consists of random extended defects, the creep expo- 
nents would change dramatically. Indeed, the latter type 
of disorder, which could occur in samples if there are, 
e.g., steps in the height of the film, allows for a variable- 
range hopping regime for creep in which the exponent 
fi — 1/3 in the two-dimensional case. Moreover, upon 
increasing the driving force, a crossover occurs to the so- 
called half-loop regime where the exponent fj, = 
An alternative explanation for the experimental results 
of Yamanouchi et alM- would be that /? 7^ so that 
the behavior for field and current-driven motion is simi- 
lar. If the pinning potential is random and extended, it 
would be possible that the current-driven experiment is 
probing the variable-range hopping regime with /i — 1/3, 
whereas the field-driven case probes the half-loop regime 
with fi — 1. This scenario would also reconcile the re- 
sults of Ref. [3^ with previous ones^i which yielded a 
critical exponent of /i ~ 0.5, as the latter could be in a 
different regime of pinning. In conclusion, further experi- 
ments are required to clarify this issue. The conjecture of 
pinning by extended defects may be experimentally ver- 
ified by increasing the driving in the current-driven case 
and checking if the exponent crosses over from ^ — 1/3 to 
/i — 1, while remaining in the creep regime. Finally, since 
the exponent /i = 1/3 occurs strictly for variable-range 
hopping in two dimensions, we note that the mapping 
presented in this paper is crucial in obtaining this result. 
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